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Two perturbation formulas have been proposed to evaluate the coupling 
between parallel optical waveguides, one involving a line integral and the 
other a surface integral. They are shown to be identical. The former 
expression is preferred because of its greater simplicity. The case of two 
parallel lossy dielectric slabs is discussed as an example. 

I. INTRODUCTION 

There has been a renewed interest during the last few years in the 
evaluation of the transverse* coupling between two parallel open wave- 
guides in connection with integrated optics circuitry 23 and long- 
distance optical communication by bundles of glass fibers. 

The coupling between two open waveguides can be obtained by 
replacing the field of one waveguide by an equivalent current and 
evaluating the perturbation caused by this current on the other 
waveguide. 4 A more direct and slightly more general (but essentially 
equivalent) derivation, based on Lorentz's reciprocity theorem, is given 
in this paper. A related result, applicable only to lossless fibers, has 
been used to evaluate the coupling between dielectric rods with circular 
cross section. 6 The perturbation formula derived in this paper involves 
an integral along a contour located between the two waveguides. A 
seemingly different perturbation formula has been recently proposed 
that involves a surface integral over the cross section. 6 The two 
formulas are shown to be in fact identical. We will not discuss in 
detail other coupling formulas such as the ones proposed in Refs. 7 
or 3. In Ref. 7, the coupling is obtained by applying the Rayleigh-Ritz 

* The word "transverse" is used here to distinguish the problem of two dielectric 
waveguides lying side by side, where the transfer of power takes place in transverse 
directions, and the axial coupling between two waveguides placed end to end, where 
the transfer of power takes place along the z axis (the later arrangement is discussed, 
for instance, in Ref. 1). 
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optimization technique to a variational expression. The formula ob- 
tained by this method involves surface integrals and is rather com- 
plicated. In Ref. 3, analytic expressions were obtained for the coupling 
between two identical rectangular fibers that agree well with numerical 
calculations based on the exact field equations. The approach, how- 
ever, is restricted to fibers with a particular geometry. 

II. GENERAL EXPRESSION OF THE COUPLING 

Let the time dependence of the sources be denoted exp ( — nt). 
Maxwell's equations are, in a source-free region with scalar permit- 
tivity e and permeability n , 

V X E = KfioH., (la) 

V X H =-/ceE. (lb) 

Any two solutions (E, H) and (E a , H Q ) of eq. (1) satisfy the relation 

V-J-O, (2a) 

where 

J = E B X H - E X H a . (2b) 

Integrating over a volume V, Lorentz reciprocity theorem is obtained 
(E XH-EXH„)-d>S = 0, (2c) 



/. 



where S denotes the surface enclosing V, and dS a vector normal to 
S pointing outward with magnitude dS. Let the medium be uniform 
along z, that is, e be independent of z. If 

(E, H) = (E„ E tj H t , H t ) exp (yz) (3) 

denotes a solution of Maxwell's equations, then 

(E+, H+) = (-#., E t , H„ -Ht) exp (-yz) (4) 

is also a solution of Maxwell's equations. The arguments x, y have 
been omitted in the above expressions, and the subscripts t stand for 
"transverse." The field (E+ H+) is the field adjoint to (E, H); it 
describes a wave propagating in an opposite direction in the same 
medium and at the same frequency. A more general definition of the 
adjoint field, applicable to nonreciprocal media, can be found in Ref. 8. 
Let us now consider two open waveguides a and b uniform along 
the 2-axis, and let S be the surface S a + £« + C a dz shown in Fig. 1. 
The field (E a , H a ) in eq. (2) is taken as the field of a trapped mode on 

218 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1974 



c„ 



C, 




Fig. 1 — Coupled dielectric waveguides. 

waveguide a in the absence of waveguide b. The dependence of (E a , H a ) 
on z is denoted exp (y a z). The field (E+ H+) is the adjoint field of a 
trapped mode of the two coupled waveguides, with an exp ( — Tz) 
dependence on z. Letting the spacing dz between <S a and S a tend to 
zero, eq. (2) becomes 

( Ta -T) f (E„ X H+ - E+ X H )-d£ a 

JS 

= - f (E„ X H + - E+ X H a ) -dCa, (5a) 
Jc a 

where &C a is a vector perpendicular to the contour C a , pointing out- 
ward. Proceeding similarly for waveguide b we obtain 

( 76 -T) f (E 6 X H+ - E + X B. h )-dS b 
Js b 

- - ( (Bi X H + - E+ X H 6 ) dC 6 . (5b) 

JCb 

Because the coupling between the two waveguides is small, we can 
assume that the field E, H at plane z = is the sum of the fields of 
the two waveguides, that is, 

E = E a + E b , (6) 

H = H a + H 6 . 
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Substituting these expressions, eqs. (6), in eqs. (5a) and (5b) we 
observe that the cross terms can be neglected on the left-hand sides 
(l.h.s.) because (E&, H 6 ) is small when (E , H„) is large, and vice 
versa. On the right-hand sides (r.h.s.), on the contrary, only the 
cross terms remain, as we can verify by applying Lorentz reciprocity 
theorem to each waveguide. Multiplying together the l.h.s. and r.h.s. 
of eq. (5a) and (5b) the desired equation for r is obtained. 

(r - 7 a) (r - y b ) = c a c b /P a P b , (6a) 

where 

Ca,6 = f (E 0l b X H 6 ta - E& X H„ , 6 ) • dC a , 6, (6b) 

JCa.b 

P a .b = f (E..6 X Hi - Ei X H Bl6 ) -dS.,6. (6c) 

JSa.b 

Because the coupling takes place only if y a ~ 76, the coupling c a 
(resp. Cb) is independent of the choice of the contour C a (resp. Cb) as 
long as it surrounds only one waveguide. By choosing the two contours 
as coincident in the region where the fields of the two trapped modes 
have a significant intensity and using eq. (4), we find that c„ is equal 
to Cb. It is shown in the appendix that our result, eq. (6), can be 
expressed in the form given in Ref. 6. The expression, eq. (6), however, 
is simpler to evaluate. 

Let us now assume that the contours C a , Cb coincide with the y 
axis and are closed at infinity where the fields vanish. The general 
expression, eq. (6), becomes 

(r - 7») (r - 76) = <?/P a P b , (7) 

= 5 / (EayHb, + E at Hby — EbyHat — Eb,H av )dy, 
P°=Jf (E«XH.)-zdxdj/, 



where 



and 



"// 



" (E 6 X H b )-zdxdy, 



where z denotes the unit vector directed along the z axis. 

Let us specialize eq. (7) to symmetrical stratified dielectric wave- 
guides such as the slabs shown in Fig. 2. The fields are assumed to 
be independent of y. For TE waves the electric field has only one 
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component E u = E. We have, from Maxwell's equations, eq. (1) 

E x = E. = H y = 0, (8a) 

H, = (w,)- l dE/dx, (8b) 

H x = - (y/Kfi )E. (8c) 

Equation (7) thus reduces to the simpler form 



(r - 7 a) (r - y b ) = (1 - Ky y *)E 2 a El/(j2 Eldx £ E\4x\ , 



(9) 



the fields E a and Eb of the uncoupled waveguides being evaluated at the 
same point between the two waveguides. 

III. COUPLING BETWEEN LOSSY DIELECTRIC SLABS 

If the waveguides are homogeneous dielectric slabs of thickness 2d 
and complex permittivity e we have 

E = exp (T710X) (10a) 

above or below the slabs and, for even modes, 

E = cosh (7 x a;)/cosh (y x d) (10b) 

within the slabs (obvious changes in the origin of the x axis were 
made). In eqs. (10a) and (10b) we have denned 



yl m k 2 - y, 


Real (y xo ) > 0, 


(11a) 


y\ = kW - y\ 




(lib) 


k 2 = v2« „ 




(lie) 



n 2 = e/e . 

The propagation constant y is known to satisfy 

7 x tanh (y x d) + y xo = 0. (12) 

(See, for instance, Ref. 6.) Substituting E from eqs. (10a) and (10b) 
in eq. (9) we obtain, using eq. (12), 

r - 7 = ±7" 1 7^>4(1 + y X od)~ l (y 2 x ~ yl)- 1 exp (-y X0 D), (13) 

where D denotes the spacing between the slabs. This expression 
coincides with the result given in Ref. 6 when appropriate changes of 
notation are made. 
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Fig. 2 — Coupled dielectric slabs. 

Let us now make a general comment. The coupling formula, eq. 
(6), rests on the existence of a divergenceless quantity, the vector J 
in eq. (2a). Coupling formulas similar to eq. (6) can be derived from 
other wave equations. For the case of the scalar parabolic wave 
equation 9 applicable to the propagation of radio waves in atmospheric 
ducts, * the vector J has components 

J, = 2k E+E « E+dE/dz - EdE+/dz, 
J t = E+VtE - EV t E+, 

where the adjoint field is 

Et{x, y, z) = E a (x, y,-z). 

The above expression for J can be obtained by analogy with the 
equivalent quantum mechanical problem. 9 

In conclusion, we have derived a simple coupling formula which is 
more general than previous similar expressions 45 because it is appli- 
cable to lossy fibers. In order to evaluate explicitly the coupling, one 
needs to know the normalized field of each waveguide, in the absence 
of the other, along some line located between the two waveguides. 
For slabs and rods with circular cross section, exact solutions are 
available. In general, however, we have to resort to numerical tech- 
niques or to measurements made at a convenient wavelength on a 
scaled version of the open waveguide. In a second part of this paper, 10 
we will apply eq. (6) to mode-selecting systems. 



* A similar equation is applicable to anisotropic fibers that have small transverse 
variation of permittivity. 8 Note that, in this approximation, a curvature of the fiber 
axis is equivalent to a constant gradient of refractive index. 
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APPENDIX 

The purpose of this appendix is to show that for lossless fibers the 
coupling formulas given in Refs. 4 and 6 coincide. 
Let (E, H) denote a field in free space 

V X E = K Mo H, 

V XH =-k€„E, 

and (E t| H&) a field in a dielectric with permittivity e(r) 

V X E b = KHollb, /,c\ 

V X H 6 =-«eEb. 

It is easy to show that these fields satisfy the relation 

f (E X Hb - Eb X H)dS = < J v (e- e )E-E b dV (16) 

in any source-free volume V bounded by S. Let now the surface S 
be the surface S b + S' b 4- C b dz shown in Fig. 1, (E 6) H 6 ) be the field 
of a trapped mode of waveguide b with an exp (jbz) dependence on z, 
and (E, H) be the adjoint field (E+, H+) of a trapped mode of wave- 
guide a, with an exp (-y a z) dependence on z. The field (E^, H^) 
satisfies eq. (14) inside the surface 5 that we have just defined. If the 
two trapped modes are degenerate, that is, if y a = 7&> the contribu- 
tions of the two surfaces S b and £» on the l.h.s. of eq. (16) cancel out. 
Therefore, letting dz tend to zero, eq. (16) becomes 



f (E+ X Hb - Eb X H+)-dCb = k f (e - e )E+-Ebd-S 6 . 

Jc b JSb 



(17) 



A similar relation can be obtained for waveguide a. Our coupling 
equation, eq. (6), can therefore be written in the form given in Ref. 6, 
except for the fact that in eq. (6) E+ and E 6 represent fields at the 
same frequency. In Ref. 6 the field E+ is defined at the opposite 
angular frequency -k, that is, E^ is replaced by E<T, where the asterisk 
denotes complex conjugation. For lossless fibers this difference is 
unimportant because E^" can be assumed real. 
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